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I. INTRODUCTION 


Collision of particles frequently take place in the accretion disks around black holes. The 
phenomenon of arbitrary high CME is a purly relativistic phenomenon and they serves as an 
excellent tool to study high energy astrophysics. Bahados, Silk and West (BSW) [it] studied the 
collision for two particles around a Kerr black hole and determined the center of mass energy 


, 3 , 


the authors further elucidated the BSW 


(CME) in the equatorial plane. Subsequently, in 
mechanism. They pointed out that the arbitrarily high CME might not be achievable in nature 
due to the astrophysical limitations i.e., the maximal spin and gravitational radiation. Lake 3,3 
demonstrated that the CME for two colliding particles is divergent at the inner horizon of a non¬ 
extremal Kerr black hole. Grib and Pavlov 3"3 showed that very large values of the scattering 
energy of particles in the centre of mass frame can be obtained for an extremal and non-extremal 
Kerr black hole. The collision in the innermost stable circular orbit for a Kerr black hole was 
discussed in 3- 33, author considered the collision for two neutral particles within the 


context of the near-horizon extremal Kerr black hole anc 


demonstrated that the CME is finite 


for any admissible value of the particle parameters. In [l]|, the authors showed that the particle 
acceleration to arbitrary high energy is one of the universal properties of an extremal Kerr black 
hole not only in astrophysics but also in more general context. An explicit expression of the CME 
for two colliding general geodesic massive and massless particles at any spacetime point around a 
Kerr black hole was obtained in [l3. They found that, in the direct collision scenario, an arbitrarily 
high CME can arise near the horizon of an extremal Kerr black hole not only at the equator but 
also on a belt centered at the equator. This belt lies between latitudes ±acos(-\/3 — 1) — ±42.94°. 

n 

In [13||, the author argued the possibility of having infinite CME in the centre of mass frame of 
collidi ng p articles is a generic feature of a Kerr black hole. 

In |14l |. the authors investigated the CME in the background of a Kerr-Newman black hole. 
They pointed out that the unlimited CME requires three conditions: (1) the collision takes place 
at the horizon of an extremal black hole, (2) one of the colliding particles has critical angular 
momentum, and (3) the spin parameter a satisfies a < 1. In 33, the author studied the 

collision of two general geodesic particles around a Kerr-Newman black hole and get the CME of 
the non-marginally and marginally bound critical particles. The collision for a freely falling neutral 
particle with a charged particle revolving in the circular orbit around a Schwarzschild black hole 


was considered in 


Ifil ]. In 17|, the authors studied the collision for two particles with different 


rest masses moving in the equatorial plane of a Kerr-Taub-NUT black hole. They demonstrated 
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that the CME depends on the spin parameter a and NUT (Newman-Unti-Tamburino) charge n. 
Exact Lense-Thirring (LT) precession and causal geodesics in the inner-most stable circular orbit 
(ISCO) in a Kerr-Taub-NUT black hole was studied in 18-20|. The CME of the collision for two 
uncharged particles falling freely from rest at inhnity in the background of a charged, rotating and 
accelerating black hole was investigated in 


In 


m 


22l |. the authors discussed the collision for two particles in the background of a stringy 


black hole. They found that the CME is arbitrarily high under two conditions: (1) the spin 
parameter o 7 ^ 0 , and ( 2 ) one of the colliding particles should have critical angular momentum. 


The collision for two particles in the background of a charged black string was discussed in 


231]. 


It was shown that the CME is arbitrarily high at the outer horizon if one of the colliding particles 
has critical charge. The particle acceleration mechanism five-dimensional compact black string, 
has been studied in Q. The. «... .he sca..e.i„. e..e... pa..ie.es i.. .he ce„.e. of 

mass frame can take arbitrarily large values not only for an extremal black string but also for 
a non-extremal black string. The CME in the absence and presence of a magnetic field around 


a Schwarzschild-like black hole was investigated in 


25(. In [2^, the authors discussed the CME 


for two colliding neutral particles at the horizon of a slowly rotating black hole in the Horava- 
Lifshitz theory of gravity and a topological Lifshitz black hole remains finite. The collision for 
test charged particles in the vicinity of the event horizon of a weakly magnetized static black hole 
with gravitomagnetic charge studied in [^. In 28], the author argued that the BSW effect exists 


for a non-rotating but charged black hole even for the simplest case of radial motion of particles 


in a Reissner-Nordstrom black hole. In 


the author gave simple and general explanation to 


the effect of unbound acceleration of particles for Reissner-Nordstrom and Kerr black holes. The 


CME of the collision for charged particles in a Bardeen black hole was studied in [30(|. In 


31], the 


authors investigated the CME near the horizon of a non-extremal Plebanski-Demianski black hole 
without NUT parameter. The CME in the background of Ayon-Beato-Garcia-Bronnikov (ABGB), 
Einstein-Maxwell-dilaton-axion (EMDA) and Bahados-Teitelboim-Zanelli (BTZ) black holes was 


investigated in 


3 |. 


A non-vacuum solution of the Einstein field equations is a Kerr-Newman-Taub-NUT (KNTN) 
black hole, which besides the spin parameter a and electric charge Q carries the NUT charge n, the 
later one plays the role of a magnetic charge. We adopt the Hamilton-Jacobi approach to study 
the dynamics of a neutral particle in the background of a KNTN black hole. We do not restrict the 
dynamics and collision to the equatorial plane alone. Instead we choose arbitrary 9 and fix 0 = ^ 
only as a special case. We discuss the detailed behavior of the CME for two neutral particles with 
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different rest masses mi and m 2 falling freely from rest at infinity in the background of a KNTN 
black hole. We derive the CME when the collision occurs at some radial coordinate r and angle 9 
close to the horizon. We show that the CME near the horizon(s) of an extremal and non-extremal 
KNTN black hole is arbitrarily high when the specific angular momentum of one of the colliding 
particles is equal to the critical angular momentum and non-vanishing spin parameter a. 

The paper is organized as follows. In Sec. II, we will discuss the equations of motion for a 
neutral particle in the background of a KNTN black hole. In Sec. Ill, we will obtain the CME of 
the collision for two neutral particles and discuss the properties. In Sec. IV, we will give a brief 
conclusion. We use the system of units c = G = 1 throughtout this paper. 


II. EQUATIONS OF MOTION IN THE BACKGROUND OF A 
KERR-NEWMAN-TAUB-NUT BLACK HOLE 


In this section, we will study the equations of motion for a neutral particle in the background of 
a KNTN black hole. Let us first give a brief review of a KNTN black hole. The KNTN black hole 
is a geometrically stationary and axisymmetric non-vacuum object, which is an important solution 
of the Einstein field equations. The KNTN black hole is determined by the following parameters 
i.e., the mass M, spin parameter a, NUT parameter n and electric charge Q. The KNTN black 
hole can be described by the metric in the Boyer-Lindquist coordinates {t,r,9,(f)) as in [33 h35| 


ds^ = —-^(A — sin^ 9)dt‘^ + -^(xA — a(S -|- ax) sin^ 9)dtd(f> + ((S -|- ay)^ sin^ 9 — x^A)d(/)^ 

2 ^ 2 ^ 2-1 


+—dr‘^ + m9^, 


( 1 ) 


where E, A and x are respectively defined by 


S = r'^ + {n + a cos 9)'^, 

A = - 2Mr - (2) 

X = a sin^ 0 — 2n cos 0. 


The KNTN metric contains the following metrics as special cases: Kerr-Taub-NUT (Q = 0), Taub- 
NUT (a = Q = 0), Kerr-Newman (n = 0), Reissner-Nordstrom (a = n = 0), Kerr (n = Q = 0) 
and Schwarzschild (a = n = Q = 0). 

The metric ([I]) becomes singular if S = 0 or A = 0, whereas E = 0 is the curvature singularity 
and A = 0 is the coordinate singularity^. Here, E = 0 implies r = 0 and cos 9 = —The 


^ The curvature invariants are given in Appendix. 






5 


horizon(s) of the KNTN black hole occur at r± = M ± where r+ and r_ 

define the outer and inner horizons, respectively, which are roots of the equation A = 0. The 
existence of the horizons require > a? + — M^, where “=” and “>” correspond to the 

extremal and non-extremal KNTN black holes, respectively. 

Now, let us discuss the equations of motion for a neutral particle of mass m in the background 
of a KNTN black hole. The motion of the particle can be determined by the Lagrangian 

^ (3) 


where the overdot denotes differentiation with respect to an affine parameter A related to the proper 
time r by r = m\. The normalization condition is ^ = At, where At = —1 for timelike 

geodesics. At = 0 for null geodesics and At = 1 for spacelike geodesics. For the massive particle, we 
have At = — 1. The 4-momentum of the particle is 


^ ^x^^ 


— 9iJ.uX , 


(4) 


which is related to the 4-velocity by 


Upi = —, (5) 

m 

dx^ 

where is the proper time for timelike geodesics. Using Eq. (BD , we can express in 

terms of the 4-momentum as x^ = g^'^Py. The Hamiltonian is given by 


n = p^x>^ -c = 


( 6 ) 


which satisfies the Hamilton equations 


dn 


dn 

dxf^ 


Moreover, the Hamilton-Jacobi equation is given by 


V = -^ = 

dX 2^ dxt^dx’^'' 


where S is the Jacobi action and 


^ = P. 
dx^^ ^ 

The Hamilton-Jacobi equation allows separation of variables in the form 


S{t, r, 9, 4>) = - £t + Ccj) + Sr{r) + Sg{6), 


(7) 

( 8 ) 

(9) 


( 10 ) 
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where £ and C are respectively the energy and angular momentum of the particle, Sr and S'e are 
arbitrary functions of r and 0, respectively. Here, and = L, where E and L are 

the specific energy and specific angular momentum of the particle defined hy E =— and L =—. 
Using these relations and Eq. ([9]), we get 

.2 


^ Pt A — sim 0 f yA — a(S + ay) sim 6 ^ 

E = -=--- u --- 

m Zj L 

r _ -Pfli _ XA - a(S + ax) sin^ (9 ^ , (S + ax)^ sin^ 6 - ^ 

1 j — — — -—- u H-—- u 

m h L 


( 11 ) 

( 12 ) 


Solving Eqs. m Qiiici fjll^p . W 6 obtctin 


u = 


X{L - xE) ^ (E + gy ) [E(E + ay) - aL] 


E sin^ 9 


AE 


^ ^ L-xE a[(E + ay)E-aL] 
E sin^ 9 AE 


(13) 

(14) 


By Eqs. ([ 8 ]) and (fTOll . we obtain 


1 /dS 0 \ 9 /,// 2 \ 4ncos0E/ \ 

—2 ( I + cos^ 0 ( (l - E^)a^ H- 5 — ) + 2ancos9(l - 2E^) -\ -— (ncos9E + L] 

\ d9 J V ^ ' sm^ 9 J sm^ 9 \ ) 

= -r^-n^-{L-aEf + ^(^{r^ + n^ + a^)E-aL^. (15) 


The left-hand side of Eq. (|15p does not depend on r while the right-hand side does not depend on 
0, hence each side must be a constant. This constant is termed as the Carter constant denoted by 
K and is a conserved quantity. Therefore 


m 


+cos^ 9( (l—E‘^)a‘^-\ - +2ancos9(l—2E‘^)-\ -^- (ncos9E+L^ = K, (16) 

^ \ d9 J \ sin 9 J sin 6* V / 


^("^) + + n^ + a^)E - aL^ --K. (17) 

Using the relations Ur = a'lid '^9 = remaining 4-velocity components are 


Eu® = ±^0, 

EvT = ±/R, 


(18) 

(19) 


with 


0 = 0 ( 0 ) = K -cos^9[{l- E'^)a^ + - 


L2 


sin^ 9 


An cos 9E f \ 

- 15 - ncos 9E + L], 

sin 9 \ / 


— 2 ancos 0 (l — 2E‘^) 


R = R{r) = {E{ax + E)-aL)^-A{K + r‘^ + n^ + {L-aE)^). 


( 20 ) 

( 21 ) 
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The ± signs are independent from each other, but one must be consistent in that choice. The +(—) 
sign corresponds to the outgoing(ingoing) geodesics. Clearly, the Carter constant K vanishes for 
the equations of motion in the equatorial plane • The radial equation of motion (fT^ can 

also be written as 




( 22 ) 


with the effective potential 

2Mr + 2n^ — sin^ 0 — + 2ancos 0 1 

2 (r2 + (n + acos2 6')2) 2(r2 + (n + a cos2 0)2)^ 

X (r^ — 2Mr — v? + + (r^ — 2Mr — (^K — a cos 6{a cos 6 + 2n)) 

T^acos 0(a^ cos^ 6 + 2acos6{r^ + 3n^ — 2ancos 9) + 4n(r^ + n^)) 

-a^ (Sn^ + + 2Mr - + 2aLE{2Mr + 2n^ - . (23) 

1 ( 7^') 

From Eqs. (fT9l) and ([22]), we conclude that Ve{{(r,9) = ^ — 1)— 2 y ;2 ■ Note that from Eqs. 

m and ()19p . for the allowed motion 0 > 0 and R>0 must be satisfied. Hence, the allowed and 
prohibited regions for the effective potential are given by 14 ff(r, 0 ) — 1 ) and 0 ) 

x(E^ — 1 ), respectively. Also, Eefr(^) 0) —)• 0 as r —oo. 

In the equatorial plane, the effective potential is given by 

T. / vr\ 2Mr + 2ri^ — 1 9 9 ^ 9 \ 

Vesir,-] = -^ + L^(r^-2Mr-n^+ Q^) 

V’ 2 y 2 (r 2 + n 2 ) 2 (r 2 + n 2 ) 2 V ^ ^ 

-a^ (3n2 + + 2Mr - Q"^) + 2aLE (2Mr + 2n2 - ^ • 


(24) 


The function E(r) can also be written in the form 


R{r) = (E^ - l)r^ + 2Mr^ + [(E^ - 1 ) 0 ^ - + 2E‘^n‘^ - - K]r^ + 2M[(L - aEf + 

+K]r + E^{3a^n‘^ + - a^Q^) + L‘^{n^ - Q^) + 2aLE{-2n‘^ + Q^) - {n^ + K) 

x(a^-n2 + Q2), (25) 


Note that coefficient of the highest power of r on the right-hand side is positive if E > 1. Only in 
this case, the motion can be unbounded (infinite). Eor E < 1, the motion is bounded (finite) i.e., 
the particle cannot reach the horizon(s) of the black hole. Eor E = 1, the motion is marginally 
bounded i.e., the motion is either finite or infinite. In this case, the particle’s motion depends on 
the black hole parameters and specific angular momentum for the allowed and prohibited regions 
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of R{r) and 0(0) but the motion can be fully analysed by R{r) or in the equatorial 

plane. The particle whose motion is bounded, unbounded and marginally bounded are respectively 
called bound, unbound and marginally bound particle. For bound and marginally bound particles, 
we have Ves{r, 9) < 0 and 14fj(r, 9) < 0, respectively. 

We need to impose the condition u* > 0 along the geodesic. This is called the “forward-in-time” 
condition which shows that the time coordinate t increases along the trajectory of the particle’s 
motion. From Eq. m, this condition reduces to 

E [(S -|- ay)^ sin^ 9 — A] > L [a(S -|- ax) sin^ 9 — yA]. (26) 


For r ^ r_|_, Eq. (|26|) implies 


L < 


E[2{n‘^ + Mr+) - Q^] 
a 


(27) 


Here, we get the upper bound of the specific angular momentum at the outer horizon of the non¬ 
extremal KNTN black hole which is called the critical angular momentum and is denoted by L_|_ 
i.e., 


L+ 


E[2(n2 + Mr+)-g2] 

a 


(28) 


Similarly, the critical angular momentum at the inner horizon of the non-extremal KNTN black 
hole is given by 

E\2(n‘^ + Mr) - 

L_ = --- - —(29) 


Eor the extremal KNTN black hole, we use r+ = M in Eq. ()28jl . which gives the critical angular 
momentum at the horizon of the extremal KNTN black hole 


L = 


E{2a^ + Q^) 

a 


(30) 


For a = 0, Eqs. ([281) . (f2^ and 1(301) become ill-defined, so we will assume a ^ 0 throughout our 
work. 


III. CENTER OF MASS ENERGY FOR TWO NEUTRAL PARTICLES 

In this section, we will study the CME of the collision for two neutral particles with different 
rest masses falling freely from rest at inhnity towards a KNTN black hole. Let us consider that 
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these particles collide at some radial coordinate r which are not restricted in the equatorial plane. 
The 4-momentum of the ith particle is given by 


Pt = (31) 

where i = 1,2 and Pj^, uf and rrii are respectively the 4-momentum, 4-velocity and rest mass (mass 
at rest at infinity) of the ith particle. The total 4-momentum of the two particles is 


PT = Pti)+Pt2y ( 32 ) 

Since the 4-momentum has zero spatial components in the center of mass frame, therefore the 
CME for the two particles is 


Pcm = -PtPtijl = + W2W(‘2))(^1^(1)m + W2W(2v)- 

Simplifying and using = — 1 in Eq. ([33]), we obtain 

Pern 


(33) 


{nil — ^- 2 )^ 


y/2mim2 


+ 1 


(34) 


/ 2mim2 

Eor the KNTN metric m, using Eqs. m, CID, (IIHD and (fl^ into Eq. (|34ll . we get the CME of 
the collision 


Er 


y/2mim2 V 2mini2 


(mi — m 2 Y ^ F{r, 9) — G{r, 9) — H (r, 9) 


I{r,9) 


(35) 


where F{r,9), G{r,9), H{r,9) and I{r,9) are given by 


F{r,9) 

G{r,9) 

R^{r) 

H{r,9) 

ei{9) 


I{r,9) 


AS sin^ 9 — {A — sin^ 9)LiL2 + ((S -|- axf sin^ 9 — x^A)EiE2 
+ (xA - a(S ax) sin^ 9) {L 1 E 2 + L 2 E 1 ), 
sin^ 9i/Ri{r)R2{r), 

[(r^ -|- TT-^ -|- a^)Ei — aLj] — A[iXj P P P (Tj — aEj)^], 


Asin2 0v^0i(0)02(0), 

Ki — 2ancos9{l — 2Ef') — cos^ “ E‘f^a? H— ^ 

4n cos 9Ei 


sin^ 9 
AS sin^ 9. 


n cos 9Ei P Li], 


(36) 


Here, Ei, Li and Ki are respectively the specific energy, specihe angular momentum and Carter 
constant of the ith particle. Clearly, the CME (I35p is invariant under the interchange of the 
quantities Li -f-)- L 2 , Ei -H- E 2 and mi o m 2 . 
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A. Near-horizon collision of particles around the non-extremal KNTN black hole 

Let us discuss the properties of the CME as the particles approach the horizons r_|_ and r_ 
of the non-extremal KNTN black hole. 


1 . Collision at the outer horizon 

The terms F{r,9) — G{r,6) — H{r,9) and I{r,9) of right-hand side of Eq. (I35p vanish at r_|_. 
Using L’Hospital’s rule and the identity — 2Mr+ — = 0, the value of the CME at 

r_|_ becomes 


Er 


y/2mim2 


(mi — m2Y drF{r, 9) — drG{r, 9) — drFl (r, 9) 


r—)-r+ 


2mim2 


drl{r,9) 


(37) 


r—>-r+ 


where 


drF{r,9)\^^^^ = 2(r_|_ — M) + (n + acos 0)^) sin^ 0 — 2(r+— M)LiL 2 


_l_ \Arj^ (r^ + + a^) sin^ 9 — 2(r_|_ — M)(a sin^ 9 — 2ncos 0 )^]Ei£'2 

—2(2n(r+ — M) cos 9 + aM sin^ 0) {L 1 E 2 + T 2 E 1 ), 


drGir,9)\ 


r—>-r+ 


sin^ 9 


_2y/ Ri{r)R 2 {r) 


R2{r)drRi{r) + Ri{r)drR2{r) 


(38) 


r— 


drRi{r)\^^^^ = 4r+ [(r^ + + a?‘)Ei - aLi]Ei - 2(r+ - M) [Ki + + (Lj - aEi )"^], 


5r-f7(r,6»)|^^^^ = 2{r+ - M)sm^ 9y^e^j9je^, 


drl{r,9)\_^^^^ = 2(r+— M)(r^ + (n + a cos 6*)^) sin^ 6*. 


After much simplification, we get the CME at the outer horizon 


Er 


2y/mim2 




(mi — m2y 

4mim2 


+ 1 + 


1 


+ 


1 


4(L+i-Li)(L+2-L2) L 
, 2\2 


[(L+i-Li)-(L+2-L2)]' 


/ f ^2 i„ 2\2 

r 2 ^ 2^ 2^2 (^1-^+2 - ^2 L+i)' + A:2(L+i - Li)^ 
\ (ri + n"' + 


+ (n + a cos 0 )^ 

+A'i(L +2 - L 2 f - a cos 0(271 + acos0) [(L+i - Li)^ + (Z /+2 - -^ 2 )^] 


1 


+ 


2 (r^ + (n + a cos 0)^) sin^ 0 

[(a sin^ 9 — 2n cos 0)^ — sin^ 0] 


cos^ 9 L 1 L 2 + 


2an cos 0 (LiL _|_2 + L 2 L^i) 


ri + 77^ + 


L+iL+2 + sin 2 0 v/ 0 i( 0 ) 02 ( 0 ) 


, (39) 


(r^ +772-1- a2)2 
r momentum foi 

The necessary condition to obtain an arbitrarily high CME is Li = L^i. 


where L+j is the critical angular momentum for the 7th particle, and can be written as L+j = 

Ai[2(772+Mr+)-Q2i 

a 
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Choosing Ei = E 2 = E, we get L+i = I /+2 = -£+ = ([39]l reduces to 


Er 


2^17111712 


r—>-r+ 


(mi — 1712 )^ 
4:77117712 


+ 1 + 


(Li - L2f 


+ 


1 


ri + (n + a cos 6Y 


4(L+- Li)(L+- L 2 ) L 

^ E^Li - L2f + K2{L+ - Lif 


+Ki{L+ - L 2 f - acos0(2n + acos0) [(L+ - Lif + (L+ - La)^] 
1 


2(r^ + (n + a cos 0)2) sin^ 9 


cos^ 6L1L2 + 2 ncos 0 £'(Li + L2) 


+ [(osin^ 9 — 271 cos 0)^ — sin^ 0] E^ + sin^ 9\/0 i(0)02(0) 


(40) 


Let us consider a marginally bound particle (S = 1) with the critical angular momentum L+. 
The conditions for the allowed region, R{r) > 0 and 0(0) > 0 give the upper and lower bounds for 
the Carter constant K given below 


< K < 

-''■min -f'-max> 


(41) 


where 


^( 1 ) = 


= 

mm 


(r + r+) 2 (r — r+) 
(r — r_) 
cos^ 0 


-r 2 - 


_|_ ^2 \ 2 


— 71 


(42) 


—-^f(a^ + + ri)^ + 4o^n^) + "^os 9 ^^2 _|_ ^^2 _|_ ^2 ^ _ 2ancos0. (43) 

sin 0 V /a sin 0 

Solving Eq. (1411) . we find that the marginally bound particle with the critical angular momentum 

reaches the outer horizon of the non-extremal KNTN black hole if the following condition is satisfied 


Ai cos^ 0 -I- i?i cos^ 0 -I- Cl cos 0 -|- Cl < 0 for any r > r+, 
of {r+r+f{r-r+) , 


(44) 


where Ai = 2a?n, Bi = + 2r^^ , Ci = 2an{a^ + 2n^ -|- 2r\) and 

Cl = —-|- (n^ -|- r^)^ -|- a^n^. 

If one chooses 0 = ^, the CME (j39p at the outer horizon of the non-extremal KNTN black hole 
reduces to 


Er 


2^17111712 


r—>-r+ 


(mi — 77l2y 
4mim2 


+ 1 + 


1 


4(L+i-Li)(L+2-T2) 


[[{L+i-Li) 


— {L+2 — -^ 2 )] + 


(LiL +2 — T 2 I/+ 1 )" 


2 M 2 -h 2 n 2 - Q2 + 2Mv'M 2 + n2 - a2 - Q2 
■{LiL+2-L2L+if 


(2M2 + 2n2 - Q2 + 2M^M‘^ + n‘^-a?-Q‘^y 


( 45 ) 
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which is indeed finite for all values of Li and L 2 except when Li or L 2 is approximately equal to 
the critical angular momentum L+j, for which the neutral particles collide with an arbitrarily high 
CME. In the case of the same specific energies, the form of the CME (145 h at r_|_ reads 

(”; - + 1 + ^ 

4mim2 4(L+ - Li)(L+ - L 2 ) V 

2M2 + 2n^ -a‘^-Q^ + 2My/M‘^ + n2 - a2 - Q2 

- -2 - E [Li - L 2 ) 

In Eiguredl we plot the effective potential of marginally bound particles for M = 1, o = 

0.8, n = 0.4, Q = 0.7211 with different specific angular momenta L = —2, — 1, 0, 1, 2.25964 
where 2.25964 is the critical angular momentum L_|_. Clearly, the effective potential is 

negative when r > r_|_, therefore the particles can reach the outer horizon. Vertical lines in the 
subplot represent the locations of the outer and inner horizons. We also plot the CME of the 
collision for Li = —2, — 1, 0, 1 and L 2 = E+. Clearly, the CME blows up at the outer horizon 
r+ = 1.00385. 




2. Collision at the inner horizon 


Similarly the terms F(r,9) — G{r,9) — H(r,9) and I{r,9) of right-hand side of Eq. (I35p also 
vanish at r_. Using L’Hospital’s rule and by simplifying the calculation, we get the CME for the 
two neutral particles at the inner horizon 


2y/mim2 


(mi — m2y 

4m 17712 


+ 1 + 


1 


+ 


1 


4(L_i - Li)(L_2 - L 2 ) L 

(^2 _j_ 


[(L_i-Li)-(L_ 2 -L 2 )]' 


(LiL_2 - L2L_iY + K2(L_i - Lif 


+ (n + a cos 9y \ (r^ + 

-|-Eri(L _2 — L 2 Y — acos9(2n + acos9) [(E-i — Li)^ + (L _2 — L 2 Y] 


+ 


2 (r^ + (n + a cos 9^) sin^ 9 

[(a sin^ 9 — 2n cos 9y — sin^ 0] 


cos^ 9 L 1 L 2 + 


2an cos 0 (LiL _2 + L 2 L- 1 ) 


(r^ -|- 77^ + 


r_ -|- 77^ + 
L_iL_ 2 + sin^ 9Y&i(d)'S>2(0) 


(47) 


This is the CME formula for the two neutral particles, where is the critical angular momentum 
at the inner horizon, which can be written as An arbitrary high CME 

can be obtained by using the condition Lj = L_i for either of the two particles. The critical 
angular momentum is same when both particles have the same specific energy and is given by 
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I N 




r 


FIG. 1: The effective potential (top figure) and center of mass energy (bottom figure) for marginally bound 
particles in the equatorial plane of the non-extremal KNTN black hole. We set M = 1, mi = m 2 = 1, 
a = 0.8, n = 0.4 and Q = 0.7211. Vertical lines identify the location of the inner and outer horizons of the 
black hole. 
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L_i = L _2 = ^- - -, while the CME ()471) reduces to 


Er 


2^17111712 


(mi — m2y 

4mim2 


+ 1 + 


+ 


1 


4(L_-Li)(L_-L2)L 

(r^ + 


(Li — L2y 


E^{Li-L 2) +K2{L_-Liy 


+ (n + a cos Oy 

+A'i(L_ — L 2 )^ — acos6(2n + acosO) [{L- — Li)^ + (L_ — L 2 )^]^ 

7 . 2n cos‘^0LiL2 + 2ncos6E{Li+L2) 
2 (r 2 + (n+ acos 6 ») 2 ) sin^(9 L 


+ [(asin^ 9 — 2n cos 0)^ — sin^ 0]+ sin^ 9^yQl{9)&2{S) 


(48) 


Let us consider a marginally bonnd particle (E = 1) with the critical angnlar momentum L_. 
The conditions for the allowed region, R{r) > 0 and 0(0) > 0 give 


< K < 


(49) 


where and iLmax are given by 


( 2 ) 


x 2 I ^2 \ 2 


All = l + r-) + r . _ 


(r - r+) 

(2) _ cos20 / 2 2 2 x2 ^, 2^ 2 


(50) 


^mL = 2 "" 2 /I (+ n‘^ + r^y + 4a^n^) + {a^ + + r^) - 2ancos0. (51) 

sm 0 V ) a sm 0 

The inequality ()49]l gives the upper and lower bounds for the Carter constant K. By Eq. (liOll . one 
can say that the marginally bound particle with the critical angular momentum reaches the inner 
horizon of the non-extremal KNTN black hole if the following condition is satisfied 


A 2 cos^ 6 + B 2 cos^ 0 + C 2 cos 0 + 112 < 0 for any r > r-, 


(52) 


where A 2 = 2a?n, B 2 = a? ^ ^ ^ — 7-2 _|_ ^^2 _j_ ^,2 _j_ 2r^ , C 2 = 2an{a? + 2 n^ -|- 2rV) and 


_ _^2 (^+^-) y ^-) I ^2^2 , m2 , ^2 x2 , ^2^2 


D 2 = —a 


(r-r+) 


+aV + (n^ + r^)^ + a^nL 


In the equatorial plane, Eq. ()47ll at the inner horizon takes the following form 


Er 


2y/mim2 


(mi — m2y 


+ 1 + 


1 


4mim2 ' 4(L_i - Li)(L_2 - L 2 ) 

l 2 , (LiL _2 — L 2 L-i)^ 


([(^- 1 -^ 1 ) 


— {L_2 — L 2 )] + 


2M2 2n2 - g2 _ 2 Mv^M 2 + n2 - a2 - Q2 

\L 1 L 2 -L 2 L 1 f 


(2M2 2n2 - Q2 _ 2M- a? - Q'^y 


( 53 ) 
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Clearly, the CME is finite for all values of Li and L 2 except when Li or L 2 is approximately equal 
to the critical angular momentum. For Ei = E 2 = E, Eq. (I53|) gives 


2y/mim2 


(mi — 1712 )^ 

4mim2 


+ 1 + 


1 


4(L_-Li)(L_-L2) 


(ii - L2f 


2M2 + 2n2 -a^-Q^- 2MJM‘^ + n 2 - a 2 - g2 
+ --E2 (Li - L 2 ) 


.(54) 


We plot the effective potential of marginally bound particles in Figure [2] for M = 2, a = 

0.6, n = 0.1, Q = 0.6 with different specific angular momenta L = —2, — 1, 0, 1, L_ where 
L_ = 2.83264. Clearly, the effective potential Ves(^r, is negative for r > r_, so the particles can 
reach the inner horizon after crossing outer horizon. The subplot shows the behaviour of 14ff^r, 
near the horizons and identify the location of the outer and inner horizons. We also plot the 
CME of the collision for Li = —2, — 1, 0, 1 and L 2 = T-. The CME is finite at the outer horizon 
and blows up at the inner horizon r_ = 1.35969. 


B. Near-horizon collision of particles around the extremal KNTN black hole 


Let us study the properties of the CME (f35|) as the particles approach the horizon of the 
extremal KNTN black hole. In the case of the extremal KNTN black hole, the NUT charge n, 
mass M, rotating parameter a and charge Q satisfies the relation v? = a‘^ + Q‘^ — M^. Using this 
relation in Eq. ()39p . we obtain 


Er 


2y/mim2 


1—>M 


(mi — 1712 )^ 

4mim2 


+ 1 + 


+ 


1 


4(Li-Li)(L2-T2)L 

(a 2 + Q 2)2 


[(Li-Li)-(L2-L2)]' 


{L 1 L 2 - L 2 L 1 Y + K 2 {Li - LiY 


+ (n + a cos 9^ V 
+Ki (L 2 - L 2 Y - a cos 0(2n + a cos 9) [(Li - hf + (L 2 - 12 ^] 


2[M‘^ + (n + a cos 0)2) sin^ 9 


cos^ 9L1L2 + 


2 ancos 9{LiL2 + L 2 L 1 ) 
2 a 2 + Q 2 


o^Lasin^^ — 2ncos0)2 — a^sin^^l „ . -—-— 

+ - (^^2 + Q2)2 -+ Sin2 0^01(0)02(0) 


• (55) 


Eq. (1551) is the CME for the two neutral particles at the horizon of the extremal KNTN black hole. 

The critical angular momentum at the horizon is given by Li = ——^-, for the iih. particle. 

The necessary condition for obtaining an arbitrarily high CME is Li = Lj for either of the two 
particles. For Ei = E 2 = E, we get the same critical angular momentum i.e., Li = L 2 = L = 
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FIG. 2: The effective potential (top figure) and center of mass energy (bottom figure) for marginally bound 
particles in the equatorial plane of the non-extremal KNTN black hole. We set M = 2, mi = m 2 = 1, 
a = 1.8, n = 0.1, and Q = 0.6. Vertical lines identify the location of the inner and outer horizons of the 
black hole. 
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E{ 2 a^+Q‘^) 


, and Eq. ([5^ takes the form 


Er 


2^17111712 


r^M 


[nil — m2y 

4mim2 


+ 1 + 


+ 


1 


A{L - Li){L - L 2 ) [ 

'{a^ + Qr ^2 


(El - E 2 )' 


^ nl2\ 2 E\Li-L2Y + K2{L-Lif 

+ [n + acosd)‘^ \ 

+Ki{L - L 2 f - acos0(2n + ocos 6 ) [{L - Lif + {L - ^ 2 )^]^ 
- 7 —^^^ -- 77 T- 71 — cos^ 6 L 1 L 2 + 2n COS OEiLi + L 2 ) 

2 (M 2 + (n + acos 0 ) 2 )sin 2 0 [ ^ i 

+ [(a sin^ 6 — 2 n cos 0 )^ — 0 ? sin^ 0 ] E‘^ + sin^ Oi/ 0 i ( 0 ) 02 ( 0 ) 


(56) 


Let us consider a marginally bound particle (E = 1) with the critical angular momentum L. 
The inequality (jlTTl reduces to 


< K < 

'^'■min — — ■''■max’ 


(57) 


where 


kYI = 2Mr + m2 - + 


— n 


(58) 


K^^'} = 

mm 


f ({a^ + 71.2 + m 2)2 + 4 a^n^'\ + ^^os V ^^2 ^2 ^ 2 -^^ _ 2 ancos 9 . ( 59 ) 

sin 0 V ) a sin 0 


4n cos 0 


Thus for the marginally bound particle with the critical angular momentum to reach the horizon 
of the extremal KNTN black hole, the following condition must be satisfied 


A 3 cos^ 0 + E 3 cos2 0 + Cs cos 0 + E 3 < 0 for any r > M, (60) 


where A 3 = 2a^n, B 3 = a^{2Mr + 3M2 + 5n2 + a2), C 3 = 2an{a‘^ + 2n2 + 2M2) and 
D 3 = —2a^rM — oyM"^ + {v? + + 0 ^ 7 ^. 

Further, if the collision occurs in the equatorial plane, the CME (j55p at the horizon of the 
extremal KNTN black hole reduces to 




2 y/mim 2 


r^M 


(mi - 7712)2 _ 

4mim2 4(Ei 


1 

Ei)(E 2 — L 2 ) 



1 — El) — ( 1/2 — E2)] 


2 


+ 


o 2 + Q 2 
(2a2 + Q2)2 


(L1L2-E2E1)' 


1 

2 


(61) 


which is indeed finite for all values of Ei and E 2 except when Li or L 2 approaches the critical 
angular momentum, for which the CME is arbitrarily high. When the electric charge Q vanishes. 





























Eq. (j61l) gives the result for the extremal Kerr-Taub-NUT black hole as obtained in Ref. 
When the specific energy of both the particles are exactly alike, then (I6ip becomes 
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Er 


2 y/mim 2 


r^M 


+ 


(mi — m2y 

4mim2 


+ 1 + 


1 


a 


4(L-Li)(L-L2) 
2 E^{Li-L2f 


(L 1 -L 2 )' 


(62) 


There must exist intervals for the spin parameter a, NUT charge n and electric charge Q to ensure 
that the marginally bound particles with the critical angular momentum L reach the horizon of 
the extremal KNTN black hole and collide at the horizon. Since the motion of the particle in 
the equatorial plane can be fully analysed by the effective potential , so with the help of 

the effective potential given in Eq. (1241) . we can determine intervals of a and n corresponding to 
different values of Q. The effective potential for the marginally bound particle with the critical 
angular momentum L is given by 


Eeff 



(r - Mf (Mr - a2 + m2 - 3^ - 

(r 2 + a2 + Q2 _ Af2)2 


(63) 


Here, the condition for the particle falling freely from rest at infinity to reach the horizon can be 
expressed as 


—^ < 0 for any r > M, 


which is equivalent to 


302 Qi 

Mr — o? + m 2-- > 0 for any r > M. 

2 2 a 2 - J - 


(64) 


(65) 


Combining with the condition 0 < n2 = a2 + Q2 _ j^‘2 gg^ M = 1, we get intervals for a and 
n for different values of Q as shown in Table HI Note that, there are two different intervals for the 
spin parameter a corresponding to the co-rotatin g an d counter-rotating orbits. For Q = 0, we get 
intervals for a and n as discussed earlier in Ref. [17]. With the increase of Q, the intervals for a 
and n become narrow. 

The maximum and minimum value of L can be obtained by the conditions 


14fr(r', ^) = 0, drV^s{r,^=Q. (66) 

Then the interval L E [Tminj T^max] can be determined from it. The intervals for the specific angular 
momentum for different values of a and Q are shown in Table Hll Note that, with the increase of 
a and Q, the interval L E [Tminj T^max] becomes wider. 
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TABLE I: The intervals for the spin parameter a and NUT charge n with different electric charge Q for the 


extremal KNTN black hole. 


Q 

a 

n 

0 

[ 1 , 72 ], [- 72 , - 1 ] 

[-i> 1 ] 

0.1 

[ 0 . 99499 , 1.40889], 

- 1.40889, - 0 . 99499 ] 


- 0.99748, 0.99748] 

0.2 

[0.97980, 1.39269], 

- 1.39269, - 0.97980] 


- 0.98974, 0.98974] 

0.3 

[ 0 . 95394 , 1.36485], 

- 1.36485, - 0 . 95394 ] 


- 0.97612, 0.97612] 

0.4 

[0.91652, 1.32389], 

- 1.32389, - 0.91652] 


- 0.95534, 0 . 95534 ] 


TABLE II: The interval L € [Amin) Amax] with different spin parameter a and electric charge Q for the 
extremal KNTN black hole. 


Q 

a = 1 

a = 1.1 

a = 1.2 

a = 1.3 

0 

[- 4.82843, 2] 


- 5.02685, 2.20171] 


- 5.2224, 2.34962] 


‘- 5.41546, 2.58869] 

0.1 


-4.83135, 2 . 01 ] 


- 5.02955, 2.06517] 


- 5.22492, 2.36131] 


‘- 5.41782, 2 . 59755 ] 

0.2 


-4.84011, 2.04] 


- 5.03766, 2.23636] 


- 5.23246, 2.39551] 


‘- 5.42488, 2.62375] 

0.3 


-4.85463, 1.82908] 


-5.05111, 2.19607] 


‘- 5.24499, 2.44978] 


- 5.4366, 2.67045] 

0.4 


-4.87481, 2.16] 


- 5.06982, 2.34748] 


‘- 5.26243, 2.53484] 


‘- 5.45292, 2.72286] 


We plot the effective potential in Figure [3] for A = — 1 and A = A in the top and 

bottom plots, respectively. Clearly, for case (III) Q = 0.3, a = 1.6, the effective potential 

is non-positive for A = — 1 but positive near the horizon r_|_ = r_ = 1 for A = A, so the 
particle cannot reach the horizon in this case for A = A. For cases (I) Q = 0.1, a = 1.2, (II) 
Q = 0.2, a = 1.3, (IV) Q = 0.4, a = 1.06 and (V) <5 = 0, a = \/2, I 4 fj^r, < 0 when r > M = 1 
for the both specific angular momenta. Hence, the particle can reach the horizon in all the four 
cases for A = — 1 and A = A. We also plot the CME of the collision in Figure 0] for Ai = — 1 and 
A 2 = A. For the case (III) Q = 0.3, a = 1.6, a does not belong to [0.95394, 1.36485], the CME 
only exists for r > 1.69668. This is because the collision for the two marginally bound particle 
with Ai = —1 and A 2 = A cannot take place at r < 1.69668. Eor the case (I), (II), (IV) and (V), 
the CME is divergent at the horizon r+ = r_ = 1. 

IV. CONCLUSION 

In this paper, we have studied the CME of the collision for two neutral particles with different 
rest masses falling freely from rest at infinity in the background of a KNTN background. Further, 
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FIG. 3: The effective potential for marginally bound particles in the equatorial plane of the extremal KNTN 
black hole for M = 1. Here, L = — 1 in the top figure and L = L in the bottom hgure. The horizon is fixed 
at r_|_ = r_ = 1. 

we have discussed the CME when the collision takes place near the horizon(s) of an extremal and 
non-extremal KNTN black hole. We have found that an arbitrarily high CME is achievable with 
following conditions: (1) the collision occurs at the horizon(s) of an extremal and non-extremal 
KNTN black hole, (2) one of the colliding particles has critical angular momentum, and (3) the 
spin parameter a ^ 0. We discovered the upper and lower bounds of the Carter constant K for a 
marginally bound particle with the critical angular momentum in an extremal and non-extremal 
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FIG. 4: The center of mass energy for marginally bound particles in the equatorial plane of the extremal 
KNTN black hole. We set M = 1 and mi = m 2 = 1. Here, Li = —1 and L 2 = L. The horizon is at 
r+ = r_ = 1. 

KNTN black hole. In the equatorial plane, we discovered that there exists intervals for the spin 
parameter a, NUT charge n and specific angular momentum L correspond to the electric charge Q 
for which not only two marginally bound particles reach the horizon of the extremal KNTN black 
hole but also the collision of these particles happens at the horizon. 
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Appendix 


The curvature invariants for KNTN metric are given by 


Xi = g^''R^v = 0, 

Z2 = = — 

X3 = = 


64Q^ 

+ (n + a cos 0)2)® 

8 (ei + 62 cos 0 + 63 cos^ 0 + 64 cos^ 0 + 65 cos^ 0 + 66 cos® 0 + 67 cos® 0 ) 

(r2 + (n + a cos 0)2)® 


(67) 

( 68 ) 

(69) 


where 


61 = 6n® — 12n®(3^ + 7n^Q^ — 90n®r^ + 12^n^Q^r^ — Mn^Q'^r'^ + 90n^r^ — 


—Qn^r^ — 6M^(n® — 15n^r^ + 15n^r^ — r®) + 12Mr(6n® — Q‘^r^ — bn^{Q^ + 4r^) 
+2n2(5QV2 + 3r^)), (70) 

62 = 4an(9n® — — ISQ^r^ — 15n^((5^ + 6r^) — 9M^(n"^ — lOn^r^ + 5r‘^) + n‘^{7Q‘^ 

+90Q^r'^ + 45r"^) + 6Mr(15n^ + lOQ^r^ + 3r^ — 10n^(Q^ + 3r^))), (71) 

63 = —2a^(—45n® + 17Q^r‘^ + 30n^(2(5^ + 9r^) + 45M^(n"^ — 6n^r^ + r^) — 3n^(7(5^ + OOQ^r^ 

+15r^) - 60Mr(6n^ + - 3n^{Q^ + 2r2))), (72) 

64 = 4a®n(30n^ + 7(5^ + 900r^ — 30M^(n^ — 3r^) — 60Mr(—3n^ + + r^) — 30n^((5^ 

+3r2)), (73) 

£5 = -a^(^-gon^-7-QOM{6n^-Q^)r +90M^{n^-r^)+ 3071"^{2Q^+ 3r^)), (74) 

66 = — 12a®n(3M^ — 3n^ + — 6Mr), (75) 

£7 = -6a®(M2-n2). (76) 


The non-trivial curvature invariants are finite at A = 0 and infinite at S = 0, hence these are 
coordinate and curvature singularities, respectively. 


[1] M. Bahados, J. Silk, S. M. West, Phys. Rev. Lett. 103, 111102 (2009). 

[2] E. Berti, V. Cardoso, L. Gualtieri, F. Pretorius, U. Sperhake, Phys. Rev. Lett. 103, 239001 (2009). 

[3] T. Jacobson, T. P. Sotiriou, Phys. Rev. Lett. 104, 021101 (2010). 

[4] K. Lake, Phys. Rev. Lett. 104, 211102 (2010). 

[5] K. Lake, Phys. Rev. Lett. 104, 259903 (2010). 

[6] A.A. Grib, Yu.V. Pavlov, Astropart. Phys. 34, 581 (2011). 



23 


[7] A.A. Grib, Yu.V. Pavlov, JETP Lett. 92, 125 (2010). 

[8] A.A. Grib, Yu.V. Pavlov, Grav. Cosmol. 17, 42 (2011). 

[9] T. Harada, M. Kimura, Phys. Rev. D 83, 024002 (2011). 

[10] A. Galajinsky, Phys. Rev. D 88, 027505 (2013). 

[11] T. Harada, M. Kimura, Class. Quantum Grav. 31, 243001 (2014). 

[12] T. Harada, M. Kimura, Phys. Rev. D 83, 084041 (2011). 

[13] O.B. Zaslavskii, Phys. Rev. D 82, 083004 (2010). 

[14] S.W. Wei, Y.X. Liu, H. Guo, C.E. Eu, Phys. Rev. D 82, 103005 (2010). 

[15] C. Liu, S. Chen, arXiv:1104.3225. 

[16] V. P. Erolov, Phys. Rev. D 85, 024020 (2012). 

[17] C. Liu, S. Chen, C. Ding, J. Jing, Phys. Lett. B 701, 285 (2011). 

[18] C. Chakraborty, Eur. Phys. J. C 74, 2759 (2014). 

[19] C. Chakraborty and P. Majumdar, Class. Quantum Grav. 31, 075006 (2014). 

[20] C. Chakraborty, arXiv:1407.6294. 

[21] 1. Hussain, Mod. Phys. Lett. A 27, 1250017 (2012). 

[22] S.W. Wei, Y.X. Liu, H.T. Li, F.W. Chen, JHEP 12, 066 (2010). 

[23] 1. Hussain, J. Phys.: Conf. Ser. 354, 012007 (2012). 

[24] A. Tursunov, M. Kolos, A. Abdujabbarov, B. Ahmedov, and Z. Stuchlik, Phys. Rev. D 88, 124001 
(2013). 

[25] M. Jamil, S. Hussain, B. Majeed, Eur. Phys. J. C 75, 24 (2015). 

[26] 1. Hussain, M. Jamil, B. Majeed, arXiv:1412.1112, DOI:10.1007/sl0773-014-2355-7. 

[27] A.A. Abdujabbarov, A.A. Tursunov, B.J. Ahmedov and A. Kuvatov, Astrophys. Space Sci. 343, 173 
(2013). 

[28] O.B. Zaslavskii, JETP Lett. 92, 571 (2010). 

[29] O.B. Zaslavskii, Class. Quantum Grav. 28, 105010 (2011). 

[30] N. Haider, Open J. Mod. Phys. 1, 1 (2014). 

[31] M. Sharif and N. Haider, Astrophys. Space Sci. 346, 111 (2013). 

[32] 1. Hussain, Mod. Phys. Lett. A 27, 1250068 (2012). 

[33] A. Grenzebach, V. Perlick, C. Lammerzahl, Phys. Rev. D 89, 124004 (2014). 

[34] D. Bini, C. Cherubini, R. T Jantzen, B. Mashhoon, Class. Quantum Grav. 20, 457 (2003). 

[35] J. G. Miller, J. Math. Phys. 14, 486 (1973). 



